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ABSTRACT 

We study the kinematics and excitation mechanisms of H2 and [Fe ll] lines in a sam- 
ple of 67 emission-line galaxies with Infrared Telescope Facility SpeX near-infrared (NIR, 
0.8-2.4/im) spectroscopy together with new photoionisation models. H2 emission lines are 
systematically narrower than narrow-line region (NLR) lines, suggesting that the two are, 
very likely, kinematically disconnected. The new models and emission-line ratios show that 
the thermal excitation plays an important role not only in active galactic nuclei but also in 
star forming galaxies. The importance of the thermal excitation in star forming galaxies may 
be associated with the presence of supernova remnants close to the region emitting H2 lines. 
This hypothesis is further supported by the similarity between the vibrational and rotational 
temperatures of H2. We confirm that the diagram involving the line ratios H2 2.121 fim/Bi-f 
and [Fe ll] 1.257yum/Pa/? is an efficient tool for separating emission-line objects according to 
their dominant type of activity. We suggest new limits to the line ratios in order to discriminate 
between the different types of nuclear activity. 
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1 INTRODUCTION 

Molecular gas emission has been detected in the inner 
few tens of parsec s of the centra l regio n of emission- 
line galaxies (e.g. Alloin & Galliano 2002; Gratadour et aL| 



this source shows that the molecular gas is distributed along the 
galaxy and i s spatially correlated with a y oung stellar population 
(30 Myr old. lStorchi-Bergmann etal1l2012l) . 



'-if 



20031: iRiffel et al. 2006a, b, 2008a; Storchi-Bergmann et al. 201 



Riffel & Storchi-Bergmannll201 lallbl) . However, the origin of these 



lines remains unclear. They appear in objects classified from star- 
burst to active galactic nucleus-dominated. In active galactic nuclei 
(AGNs), unification models predict the presence of a dusty molecu- 
lar torus that obscures the central source from some viewing angles, 
thus creating apparently different AGN classes (e.g. Seyfert 1 and 
Seyfert 2 galaxies). In order to block the radiation coming from 
the central source, the torus must have a very large optical depth 
(tens of magnitudes), thus serving as a natura l reservoir of molecu- 
lar gas. For example. iGratadour et al.l (I2003L see their Fig. 4) pro- 
pose that the torus is the source of the H2 emission-lines observed 
in NGC 1068. However, Integral Field Unit (IFU) spectroscopy of 
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Observational evidence has confirmed that circumnu- 
clear/nuclear starbursts can coex is t in objects harbouring 
A GNs (e.g. iMizutani et al.l 



Sturm et al 



coex is t in objects harbouring 
I1994I; llmanishi & Dudley! |200(I; 



_ .."2002l;IContini e t al."2002'; 'Rod riguez- Ardila & Viegasl 

20031; llmanishi & Wadal 12004 ; Storchi-Bergmann et all 120051: 
Shi et alj I2OO6I) . However, in the case of low-luminosity 



AGNs, no strong starburst evidence is found, but significant 
fracti o ns of intermediate-age s t ars are detected (ilv ano v et al] 
2OOOI; Icid Femandes et al.l I2OO4I ; iGonzalez Delgado et al] |2004 
Cid Femandes et alj|2005ll2010l and references therein). Thus, if 
star formation is taking place or occurred recently, the presence 
of molecular clouds associated with the star forming regions 
should leave spectral signatures that may contaminate the AGN 
spectrum. Therefore, whether the H2 lines observed in AGNs 
are directly associated with either the molecular dusty torus or 
the circumnuclear gas remains an open question. Alternatively, 
the molecular gas may be directly associated with the AGN, but 
distributed within the circumnuclear region. This later scenario is 
supported by the fact that the H2 gas is most probably arranged 
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in a disc surrounding the nucleus d Reunanen et"ai]|2002l l2003l : 
iNeumaver et alj l2007t) . In this picture, shocks by radio jets and 
X-ray heat ing are plausible excitation mechanisms for the H2 
molec ule jKnop et al] Il99€l : iReunanen et alj I2OO2I : lOors et alj 
l2012h . 

As shown by [Rodriguez- Ardila et al. ] ( |2004 l2005h and 
iRiffel et alj ( l2006ij) in the study of 51 objects, mostly AGNs, 
molecular hydrogen is detected in almost all Seyfert 2 (Sy 2) and 
in about 80% of the Seyfert 1 (Sy 1) sou rces. Interestingly, the four 
starburst galaxies in iRiffel et al.l ( l2006al) also show the presence of 
molecular hydrogen emission. However, these previous works have 
not included a statistically significant number of Low-ionisation 
nuclear emission-line region (LINER) and Star Forming Galaxies 
(SFGs). In fact, up to now, only a few investigations have included 
these type of sources (e.g. iLarkin et al. 1998). These classes de- 
serve separate analyses, as it is expected that the active nucleus has 
very little influence on their emission-line spectrum. 

The H2 infrared emission lines may be produced by the ra- 
diative de-excitation following collisional excitation (thermal pro- 
cess) or UV (6-18 eV) photon absorption (non-thermal process). 
H2 may also be formed in an excited state, after which de- 
excitation also leads to emission of IR lines. This mechanism is, 
however, not ver y important in a hot or stro ngly UV/X-ray ir- 
radiated gas (e.g. lAleman & Gruenwaldl 1201 ll) as AGNs. H2 IR 
emis sion can be significant in regions illuminated by UV p ho- 
tons telack & van Dishoecklll987l:]stemberg & Dalgamolll989l) or 



X-ravs ( lMalonevetal.lll996l:lLePD&McCravlll983l) , as well as 
in shocks i Hollenbach & McKej 198sl ). In each case, a different 
emission-line spectrum is expected and, therefore, the relative in- 
tensities among the H2 emission-lines may be used to d iscrimi- 
nate between the powering mechanisms (e.g. lMourilll994l) , keep- 
ing in mind that multiple emitting regions and multip le mech- 
anisms may simultaneously be present ( Rodriguez- A rdila et al 



200' 



)e pri 



2004 2005 ; Riffel et al. 2006bl l2008al: IStorchi-Bergmann et al 



Riffel & Storchi-Bergmannll2011 



iRodriguez-Ardila et al. I ( |2004 l2005h carried out 



a pio- 
neering study of the H2 excitation mechanisms in a sample of 
Seyfert galaxies. The i nnovation of these i nvestigations com- 
pared to previous ones ijVeilleux et al. I ll997t lLarkinetal.|[T99^ : 
IReunanen et al.ll2002l |2003|) is that their integrated spectra cover 
the inner 300 pc in most of the sample, minimizing the host 
galaxy contamination and maximizing the number of diagnostic 
lines observed. In addition, [Fell] and H2 emission-lines are 
usually suggested to form in the same region, and the evidence 
accumulated shows that forbidden iron emission in AGNs can 
have different sour ces, but very l i kely, all directly related to 
the central engine jForbes & Ward Il993l: iGoodrich et all Il994 
Simijson etal] Il99d~ Alonso-Herrero et ahf 1997t Mouri et al I 



2OOOI: 



Rodriguez- Ardila et al.l 
Riffel etal.' 2006t/, '2008;/; 'Storchi-Bergman n et alT 



12005 



Davies et al. 



2004 



2005; 



200s 



Riffel & Storchi-Bergmann 2011b; Dors et al. 20121 and ref- 
erences therein). Thus, it is interesting to study the excitation 
mechanisms leading to H2, together with the processes responsible 
for the [Fe II] emission-lines. 

With the above in mind, we present a study of the excitation 
mechanisms leading to the emission of H2 and [Fe II] for a sample 
of nearby LINERs and SFGs. We present 16 new spectra acquired 
with the sam e ins trumentation employ ed bv lRodriguez-Ardila et al.l 
( I2OO4I2OO5I) and IRiffel et"al] ( l2006j) , with simultaneous observa- 
tions of the J H K-hw\As, eliminating differences in aperture and 
seeing across the bands. The samples combined correspond to an 
homogeneous data set with ~ 70 sources. To date this is the largest 



number of galaxies used to study the nuclear excitation mecha- 
nisms of H2 and [Fe II]. 

This paper is structured as follows. § |2| presents the obser- 
vations and data reduction steps. § [3] discusses the kinematics of 
the H2 and [Fe II] gas inferred from their line profiles. In § |4] 
the excitation mechanisms of the H2 IR lines in AGN are dis- 
cussed using new photoionisation models for AGNs and Starbusts. 
The calculated emission is compared to the observed H2 line ra- 
tios in a 1-0 S(2)/l-0 S(0) vs 2-1 S(l)/l-0 S(I) diagnostic dia- 
gram. § |5] analyses the line ratios H2 2. 121/im/Br7 and [Fell] 
1.257/im/Pa/3 for different emission-line objects and the apparent 
correlation between these ratios over a large range of values. It also 
discusses the difference in these ratios for different nuclear activity 
type. Concluding remarks are given in § |6| A Hubble constant of 
75 km s^^ Mpc^^ will be used throughout this work in order to be 
consistent with our previous studies. 



2 OBSERVATIONS AND DATA REDUCTION 

Near-infrared (NIR) spectra in the range 0.8-2.4 ^m were ob- 
tain ed on October 4,6, and 7 in 2010 with the SpeX spectrome- 
ter ( iRavner et aT]|2003h attached to the NASA 3 m Infrared Tele- 
scope Facility (IRTF). The detector is a 1 024 x 1024 ALADDIN 3 
InSb aiTay with a spatial scale of 0. 15"/pixel. Simultaneous wave- 
length coverage was obtained by means of prism cross-dispersers. 
A 0.8" X 15" slit was used during the observations, giving a spec- 
tral resolution of 360 km s^^ . Both the the arc lamp spectra and the 
night sky spectra are consistent with this value. The seeing varied 
between 0.4"-0.7" over the different nights. 

Due to the extended nature of the sources, the observations 
were done by nodding in an Object-Sky-Object pattern with typical 
individual integration times of 120 s and total on-source integration 
times between 18 and 58 minutes. During the observations, AO V 
stars were observed near each target to provide telluric standards at 
similar air masses. They were also used to flux calibrate the sample. 
Table[T]shows the observation log. The galaxies are listed in order 
of right ascension, and the number of exposures refers to on source 
integrations. 

Spectral extraction and wavelength calibration were per- 
formed using the SPEXTOOL, a software developed and provi ded by 
the SpeX team for the IRTF community ijCushing et al."200 3tEi fol- 
lowing the same procedures as described by .Rodrf guez-Ar dila et al.l 
( I2OO4 ). Column 10 of Table [T]lists the radius of the integrated re- 
gion, with centre at the peak of the continuum light distribution for 
every object of the sample. No effort was made to extract spectra 
at positions different from the nuclear region, even though some 
objects show evidence of extended emission. Telluric absorption 
correction and flux calibration were applied to the individual I-D 
spectra by means of the IDL routine xtellcor jVacca et alj2003h . 

The 1-D wavelength and flux calibrated spectra were then cor- 
rected for redshift, determined from the average 2 measured from 
the position of [Slll] 0.953^m, Pa(5, He I 1.083/im, Pa/3 and Br7. 
Final reduced spectra, in the spectral regions of interest in this 
work, are plotted in Figures [T]and|2] 



^ SPEXTOOL is available from the IRTF web site at 
|http://irtf.ifa.hawaii.edu/FaciLity/spex/spex.html| 
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Table 1. Observation log and basic properties of the sample. 



Source 


Activity 


Activity 


Obs. Date Exp. Time Airmass 


RA 


DEC 


z 


Pos. Angle radius'' 






Reference 


(s) 








(deg) (pc) 



NGC23 


SFG 


1 


2010 10 07 


29 X 120 


1.04 


00h09m53.4s 


+25d55m26s 


0.0157202 


330 


609 


NGC 520 


SPG 


2 


2010 10 04 


16 X 120 


1.04 


01h24m35.1s 


+03d47m33s 


0.0080367 


300 


311 


NGC 660 


LINER/Hll 


2,3,4 


2010 10 06 


24 X 120 


1.01 


01h43m02.4s 


+13d38m42s 


0.0029152 


33 


107 


NGC 1055 


LINER/Hll 


2,3,4 


2010 10 04 


16 X 120 


1.07 


02h41m45.2s 


+00d26m35s 


0.0036267 


285 


210 


NGC 1134 


SFG 


5 


2010 10 04 


16 X 120 


1.11 


02h53m41.3s 


+13d00m51s 


0.0129803 





503 


NGC 1204 


LINER 


6 


2010 10 07 


16 X 120 


1.23 


03h04m39.9s 


-12d20m29s 


0.0154058 


66 


597 


NGC 1222 


SFG 


7 


2010 10 06 


24 X 120 


1.13 


03h08m56.7s 


-02d57ml9s 


0.0082097 


315 


270 


NGC 1266 


LINER 


7 


2010 10 07 


18 X 120 


1.09 


03hl6m00.7s 


-02d25m38s 


0.0077032 





298 


UGC 2982 


SFG 


8 


2010 10 04 


9 X 120 


1.11 


04hl2m22.4s 


+05d32m51s 


0.0177955 


295 


689 


NGC 1797 


SFG 


1 


2010 10 07 


16 X 120 


1.23 


05h07m44.9s 


-08d01m09s 


0.0154111 


66 


597 


NGC 6814 


Sy 1 


7 


2010 10 07 


16 X 120 


1.17 


19h42m40.6s 


-10dl9m25s 


0.0056730 





220 


NGC 6835 


7 




2010 10 06 


22 X 120 


1.21 


19h54m32.9s 


-12d34m03s 


0.0057248 


70 


166 


UGC 12150 


LINER/Hll 


9 


2010 10 04 


15 X 120 


1.08 


22h41ml2.2s 


+34dl4m57s 


0.0214590 


37 


748 


NGC 7465 


LINER/Sy 2 


10 


2010 10 06 


12 X 120 


1.03 


23h02m01.0s 


+15d57m53s 


0.0066328 


340 


257 


NGC 7591 


LINER 


7 


2010 10 07 


16 X 120 


1.03 


23hl8ml6.3s 


+06d35m09s 


0.0165841 





642 


NGC 7678 


SFG 


11 


2010 10 04 


16 X 120 


1.01 


23h28m27.9s 


+22d25ml6s 


0.0120136 


90 


419 



Table Notes: SF G: Star Forming Gal axie s (Star burst orHlI gal axies). LINER/HII w ere assumed as pure LI NE Rs in the text. References - 1 : | Balzanol fl983l): 2: 
IHo et alj fl997bl): 3:IHo et alj <1997al): 4:lFilho et al. (2004); 5:ICondon et al.H2002h: 6: ISturm et aljh006l) : 7: |Pareira-Santaella et al.ll2010l) : 8: ISchmitt et alj 
<200d) : 9: Tveilleux et alJil995h : 10: lFerruit et al.l i l2000i) : ll: lGoncalves et aljh997bl) . 



3 ON THE MOLECULAR GAS KINEMATICS 

The main goal of this section is to discuss how the velocity width 
of the H2 2.1213/im emission-line compares to that of [Fell] and 
[Slll] in different objects. This provides important constraints to the 
location of the molecular gas, and allows a comparison between 
this region with the different types of nuclear activity. It is worth 
mentioning that the nuclear activity used here (Table [T) comes pri- 
marily from optical studies, and was taken from the NASA Extra- 
galactic Database (NED). However, NIR emission-line ratios may 
reveal the presence of "hidden" AGNs (see Sec.|5)- 

Previous results, mainly on Seyfe rt galaxies, are con trover- 
sial. An imaging study of H2 emission jOuillen et al.|| 19991) shows 
that the H2 is distributed on scales of a few hundred parsecs from 
the nucleus. In addition, the fact that the H2 emission was found 
to coincide with [Olll] and Hq-i-[Nii] for some sources led the 
authors to suggest that the molecular gas ma y follow the n arrow 
line region (NLR) gas distribution. Schinner er et al.l (Il999h sug- 
gest that the molecular gas in NGC 1068 and NGC 3727 origi- 
nates from a warped disc with a radius sma ller than ^ 75 pc. These 
latter r esults are supported by those of iRodriguez-Ardila et al.l 
( |2004 l2005h . who found that the H2 2.1213/im emission is 
spectroscopically unresolved or with FWHM systematically nar- 
rower than that of the NLR forbidden lines in Seyfert galax- 
ies. This was interpreted in terms of a kinetic disconnection 
between the molecular and the NLR gas. Further support to 
this hypothesis comes from Integral Field Unit (IFU) studies of 
AGNs, which suggest that the H2 is in the galaxy plane, from 
the centre up to the field limit (~ 500 pc), while the ionised 
gas is observed up to h igh latitudes 1 Riffel et al ]|2006bl, l2008d; 
Storchi-Bergmann et all [20 09. 20 10l: iRiffel & Storch i-Bergmann 



201 1 al ibi: iDavies et alJbOOgljMuller Sanchez et al.ll2'o09l) . 



In order to study the II2 and [Fe II] lines in LINERs and SFGs, 
and to complete our analysis of these transitions in emission-line 
objects, we have followed the methodology of our previous stud- 
ies. We assume that the widths of the forbidden lines (or NLR gas 
in the case of AGNs) reflect the large scale motions of the emitting 



clouds in the gravitational potential of the central mass. As a con- 
sequence, similar FWHM of molecular and atomic forbidden lines 
would indicate that these species are co-spatiaQ. 

TablelUhsts the intrinsic FWHM of [S III] 0.9531/im, [Fell] 
1.2567/im and H2 1.957, 2.0332, and 2.121/im. To obtain the in- 
trinsic FWHM, the instrumental width was subtracted from the 
observed line width in quadrature. The errors in FWHM (~ 
30 kms^^) are dominated by the uncertainty in the continuum 
placement. To ensure that the intrinsic FWHM is equal or larger 
than the instrumental one (360 kms'^), a line was considered to 
be spectroscopically resolved if its measured FWHM was larger 
than 500 kms"\ Thus, lines with FWHM equal to 360 kms"^ 
in Table |2| may have measured values in the range 360 kms"^ ^ 
FWHMsC 500kms"\ 

Figure|3]compares the observed line profiles of [S III] 9531 A 
(dotted line), [Fell] 1.257/im (dashed line), and H2 2.122/im 
(full line) for a sub-sample of our sources (those exhibiting all 
the three lines). Clearly, [S III] and [Fell] tend to be broader 
than H2, following the same trend observed in Seyfert galaxies 
jRodn'guez-Ardila et al.l2004 l2005h . 

To better understand the kinematics of the gas emitting the 
different lines in the sources of different classification, we show 
in Figure |4| a fractional histogram of the FWHM values of the 
three lines used in Figure [3]for the data presented in Tabl e |2]plu s 
literature data taken from iRodriguez-Ardila et al. ] (I2OO4I l2005h . 
This plot shows that the molecular gas follows different kine- 
matics than that followed by the atomic forbidden emission-line 
gas not only in Seyferts, but also in LINERs and SFGs. Previ- 
ous work supports this disconnection in Seyfert galaxies, and sug- 
gests that the gas is arranged in a disc suiTounding the nuclear re- 



^ It is woith mentioning that, only because two lines have similar FWHMs 
does not prove that they are co-spatial. For instance, they could originate 
approximately in the same region but from different gas layers with different 
physical conditions. Even when lines originate froin exactly the same gas, 
different line profiles (or FWHM) may result from radiative transfer effects. 
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Figure 1. SpeX final reduced spectra, in the Earth's frame of reference, centred near Pa/3 (left panel), the H-band (1.60/im, middle panel), and Br7 (right 
panel). The observed flux is in units of 10~^^ erg cm~^ s~^ A^^. The identified [Fe 11] (dashed lines) and H2 lines (dotted lines) are marked in the spectra. 



gion ( iReunanen et alj2002l . l2003l) . Furt her support for this hypoth- 
esis comes from Gemini-IFU studies jRiffel et alj|2006bl 12008 j; 
Storchi-Bergmann et all I2OO9I I2OI0I : iRiffel & Storchi-BergmannI 
2011al lbh. who suggest that the H2 is located in the galaxy plane 
and is distributed along the whole field, while the ionised gas is 
observed up to high latitudes from the galaxy plane. 

In summary, the observational evidence presented here con- 
firms that H2 is common within the inner few hundred parsecs of 
emission-line galaxies, regardless of type. The molecular gas fol- 
lows different kinematics than that of the ionised gas, suggesting 
that the two emissions are not co-spatial. The H2 lines are naiTower 



than the forbidden ones, in particular [S III]. A possible exp lana- 
tion involves atomic gas outflow. Indeed, iRiffel et alj l2012h have 
shown very recently that the velocity dispersion of gas restricted to 
the plane (H2 gas) does not correlate with that of bulge stars. On 
the other hand, the velocity dispersion of gas extending to higher 
latitudes (such as the [Fe II] emitting gas) is similar to that of the 
galaxy bulge stars. Higher velocity dispersions of [Fe II] relative to 
stars are probably due to an extra heating provided by a nuclear 
AGN outflow. In such a scenario, the atomic gas could be farther 
out from the supermasive black hole (SMBH) than the molecular 
gas. Another possibility worth mentioning is the atomic gas be- 
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Figure 2. Same as Fig.[T]for the remaining objects. 



ing closer to the SMBH than the molecular gas and, thus, being 
more affected by the gravitational pull of the putative SMBH in 
AGNs. Such a possibility is related to the fact that [S III] is a higher 
ionisation line and, thus, its bulk should be formed more inwards 
in the NLR than [Fell] and H2. Thus, the unresolved FWHM of 
H2 in almost all objects implies that the molecular gas is proba- 
bly not gravitationally bound to the SMBH, but to the collective 

i yavitational potential o f the g alaxy. According to literature results 
Reunanen etal]|2002l . |2003|) . it may be arranged in a disc-like 
structure. Further support of this hypothesis comes from the fact 
that, compared to SFGs, AGNs tend to have significantly broader 
[S III] than [Fell] and H2 (Fig.©. 



4 EXCITATION MECHANISMS OF THE H2 LINES 

As mentioned in the Introduction, H2 line ratios are used to dis- 
criminate among thedominant line excitation mechanisms. For 
example, iMouri ( 1 19941) suggested the use of the ratios 1-0 S(2) 
2.033/im/l-O S(0) 2.223Atm versus 2-1 S(l) 2.247^im/l-0 S(l) 
2.121/xm to separate the mechanism exciting the H2 emission. Fig- 
ure |5lshowsJhi^_2l£t_^o^_te objec in Table [T]plus AGNs 
from TRodrfguez-Ardila et aL I(l2004l l2005h and blue compact dwarf 
galaxies (BCGs) from llzotov & ThuanI feOllh . It can be seen from 
Fig. [5] that objects with different nuclear activity populate differ- 
ent regions of the diagram. Syls and SFGs tend to scatter over 
the plotted region, although Syls tend to have lower values of 2- 
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Table 2. FWHM (in kms~^) corrected for instrumental resolution, this 
correction was applied for lines with FWHM > 500 kms~^. Lines with 
smaller values were reported as unresolved (FWHM=360 kms~^). 



Source 


[S III] 


[Fe II] 


H2 


H2 


H2 




9530 


1.2570 


1.9570 


2.0332 


2.1213 




A 


/xm 


/xm 


lim 


/jm 


NGC23 


435 


563 


492 


360 


360 


NGC 520 


_ 


_ 


360 


360 


360 


NGC 660 


612 


402 


390 


360 


360 


NGC 1055 












NGC 1134 












NGC 1204 


562 


301 


360 


360 


360 


NGC 1222 


500 


447 


601 


360 


360 


NGC 1266 




632 


360 


360 


360 


UGC 2982 






648 


360 


360 


NGC 1797 


573 


285 


360 


360 


360 


NGC 6814 


717 


583 


368 


360 


360 


NGC 6835 




499 


876 


360 


360 


UGC 12150 


611 


508 


527 


360 


360 


NGC 7465 


557 


492 


360 


390 


360 


NGC 7591 


669 


367 


666 


363 


360 


NGC 7678 


746 
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Table Notes: The uncertainty associated with the FWHM is ' 
(for details see text). 



^30 km s 




Figure 3. Observed line profiles of [S III] 953lA (dotted blue line), [Fell] 
\.251^m (dashed green line) andH2 2.122/xm (full line). SFGs are labeled 
in blue color 



1 S(l)/l-0 S(l) than SFGs, while Sy2 and LINERs are concentrate 
in a smaller region around the coordinates (0.1,1.5). Note also that 
two BCGs are on the lower right of the plot. 

Figure |5] also shows results of model calculations found 
in the literature. The UV-e xcited low-density gas models of 
iBlack & van Dishoeckl ( 1 19871) occupy the region indicated by the 
orange fill ed circle on the right, while the UV-excited high- density 
models of IStemberg & Dalgamol ( Il989h occupy the region indi- 
ca ted by the green box on the left. IVIodels of X-ray irradiated gas 
0f lLepD&lV[cCravl ( ll983h a re indicated by purple stars. The cyan 
star is a model of shock by iKwanI {191% . llMalonev et~al] jl99d) 
showed that for X-ray-irradiated molecular gas, the 2-1 S(l)/1- 
S(I) ratio is <0.3 (~0.1 in regions where the line 1-0 S(l) is more 
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Figure 4. Histogram s howing the FWHM distri bution for LINERs, SFGs 
and Sys (the latter from iRodriguez-Ardila et al...2 004, 2005 ). The measu re- 
ments for Seyferts were taken from lRodriguez-^dila et al . H2004l2005t) . 



intense). [Dorset al.l ( f2012l) calculated the ratios with the photoion- 
isation code CLOUDY for AGN models. They found that heating 
by X-rays produced by active nuclei is a common an d very impor- 
tant mechanism of H2 excitation. Shock models bv lKwanl(T977h 
give 2-1 S(1)/1-0S(1) = 0.09 2 and 1-0 S(2)/l-0 S(0) = 0.87, while 
iHollenbach & IVlcKeel [l989l) find ratios in the range 0.3-0.5 for the 



2-1 S(l)/l-0 S(l)rati(j 

The observed ratios are not reproduced by any of pure low- 
density UV, high-density UV, X-rays or shock models only. Differ- 
ent processes may be more important in different objects. IVIore- 
over, a mixture of these processes are likely to be exciting the 
molecule in both AGNs and SBs. Here we will focus our discussion 
in the excitation due to a ra diation source (i.e., no sh oc k exci tation). 

The models of IStemberg & Dalgam 

3 ( Il989l) and 

iBlack & van Dishoeckl jl987h only consid er UV radiation (E < 
13.6 eV). On the other hand, the models of llVIalonev e t al.' fl99^ 
only c onsider X-ray radiation (E > 100 eV). iLeix) & IVIcCravl 
( 1 19831) present very complete models and take a wi der energy 
range into account, but they only calculated two models. lOors et al] 
(2012J obtained models for a single value of gas density {uh = 
10^ cm~'^) and their calculations stop when the temperature falls 
to T = 1000 K, missing an important part of the H2 emission, 
where the X-rays may be of great importance (s ee the discussion 
of Fig.[8]below and lAleman & Gruenwaldll201 ih . 

Here we present new calculations of the H2 IR emission lines 
produced with a photoionisation numerical code. These calcula- 
tions are an improvement over the previous mentioned calcula- 
tions, since they consider the spectrum from the IR up to X-rays 
(10~^ keV ^ E < I keV), they run from the inner and more ion- 
ized zones to the outer and more neutral zones, where T decreases 
to lOOK. The processing of the radiation in the inner gas shells is 
done in a self-consistent way and the microphysics of H2 is in- 
cluded in detail. The models are described below (Sect. |4Ti . In 
Sect. 14.21 we compare model and observed H2 line ratios, with the 
aim of studying the excitation mechanism of these lines. 



This ratio is not indicated on Fig. |5]because the authors did not provide 
the con-esponding 1-0 S(2)/l-0 S(0) ratios. 
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Figure 5. Objects with different nuclear activity populate different regions of tlie H2 1-0 S(2) 2.033/im/l-O S(0) 2.223^m versus 2-1 S(l) 2.247/^m/l-O S(l) 
2.121/^m diagnostic diagram. Different types of nuclear activity are shown in each plot as follows: Syl, QSO, and NLSyl (black open squares); Sy2 (red filled 
squares); LINERs (green asterisks); SFGs (blue filled circles) and BCGs (purple open circles). The orange dashed curve represents the ratios for an isothermal 
and unifonn density gas distribution; each dot represents the tempera tures from 1000 to 6000 K, in steps of 1000 K, from left to light. The orange circle on 
the right covers the locus of the non-thermal UV excitat i on mo dels of lBlack & van Dishoeckl fl987l) . whi le the green box (on the left) covers the locus of the 
thermal UV exc i tation models oflStemberg &Dalgamol ( ll989l) . Purple stars are thermal X-ray models of lLepp & McCravl(ll983h and cyan star is a model of 
shock bv lKwanI h977h . lDors et alj )2012h are indicated by the two pink polygons. 



4.1 Models 

We calculated the intensity of H2 emission lines of AGNs 
using the one-dimensiona l photoionisation code AANGABA 
dOmenwald & Viegas|[T992l) . The H2 microphysics is included in 
detail in the code. To calculate the intensity of H2 IR lines, the den- 
sity and the population of the H2 electronic ground state rovibra- 
tional levels that originate the emission must be known. The density 
of H2 is calculated assuming the chemical and ionisation equilib- 
rium between this molecule and the other H-bearing species (H'', 
H+, H", , and K^; see Aleman & Gmenwald 2004). The code 
includes over forty different reactions of formation and destruction 
of these species. The population of the H2 rovibrational levels of 
the three lowest electronic bound states is calculated by assuming 
statistical equilibrium (i.e., the total population rate of a level is 
equal to the total depopulation rate). The populations of the elec- 
tronic ground level takes into account the radiative and collisional 



excitation and de-excitation mechanisms, as well as the possibil- 
ity that H2 is produced or destroyed in any given level. For upper 
electronic levels, only radiative electronic transitions between each 
upper state and the ground state are included, since this must be the 
dominant mechanism. 

The population of excited H2 rovibrational levels of the elec- 
tronic ground state occurs mainly by electric dipole transitions to 
upper electronic states, with the absorption of a UV photon, fol- 
lowed by the subsequent decay to the ground state (with the emis- 
sion of an UV photon), usually in a excited rovibrational level. 
This mechanism of H2 excitation is called UV pumping. The de- 
cay to lower levels, through quadrupole transitions, produces fluo- 
rescence. Collisions with the dominant species of the gas may also 
be an important route of H2 excitation. We included collisions of 
H2 with the main components of the gas, i.e., H, H+, He, H2, and 
electrons. 
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Table 3. Model elemental abundance relative to H 

Element Abundance Element Abundance 



to the density of H atoms of the cloud. The ionisation parameter ([/) 
is then 



H 


1.0 


Mg 


3.8x10-^5 


He 


9.8x10-2 


Si 


3.5x10-5 


C 


3.6x10-* 


S 


1.6x10-5 


N 


1.1x10-" 


CI 


3.2x10-^ 


O 


S.SxlO-* 


Ar 




Ne 


1.2x10-" 


Fe 


4.7x10-5 



The H2 molecule is also included in the gas temperature calcu- 
lation, which assumes thermal equilibrium (the total input of energy 
in the gas is balanced by the total loss of energy per unit time and 
volume). The relevant mechanisms of gain and loss of energy by 
the gas due to atomic species, dust, and H2 are taken into account. 
The gas heating mechanisms are photoionisation of atoms, atomic 
ions, and H2 by the primary and diffuse radiation, H2 photodis- 
sociation (direct and two steps), H2 formation on grain surfaces, 
by associative detachment, and by charge exchange with H, H2 
collisional de-excitation, and photoelectric effect on dust surfaces. 
The gas cooling mechanisms are emission of collisionally excited 
lines; radiative and dielectronic atomic recombination, thermal col- 
lisional atomic ionisation, free-free emission, collisional excitation 
of H2 destruction of H2 by charge exchange with H'*' , H2 colli- 
sional dissociation, and collision of gas-phase particles with dust 
grains. More details on th e calculation of the H2 level population 
are found in |Aleman & G ruenwald (2004, 201 1). 

The emissivity of a line produced by de-excitation from a 
given level is the product of the population of that level by the Ein- 
stein coefficient of the transition and by the energy of the emitted 
photon. The emissivities of the H2 IR lines are calculated for each 
point of the nebula. The intensities are calculated by integrating 
the emissivity over the nebular volume. The gas is assumed to be 
optically thin for the H2 lines. 

We assume a spherical distribution, uniform density and com- 
position for gas and dust, and gas d ensities jnu) fro m 10'' to 10''' 
cm^'*, which are typical for AGNs ([Osterbrockll 19891). Solar abun - 



dance is assumed for the elements in a ll models jOsterbrockll98S ), 
with values (Table O obtained from lOrevesse & Anders jl 989 ). 
The chemical composition affects the processing of the radiation 
inside the cloud and, consequently, the ionisation stiucture of the 
cl oud, but this does not significantly affect the H2 ratios, as shown 
bv lDorsetaLl ( l2012h . 

Dust, in our models, is assumed to be compose d by spherical 
raphi te grains with op t ical p roperties taken from lOraine & Led 
19841) : E aor & Drainj l ll993h . Models were calculated for two 
different grain radii, 10"^ and ICT^ /im, typical i nterstellar val- 
ues according to i Li 20071) a nd lMathis et al.l ( Il977h . As shown in 
lAleman & Gruenwaldl ( l20()4l) . the dust particle size is more impor- 
tant for the H2 density than the grain type. The dust-to-gas mass 
ratio (Md/Mg) is assumed constant in all models. Models were cal- 
culated for Md/Mg = 10"^ and 10"^ , average values for planetary 
nebulae and the interstellar mediu m jTielensI [20051 : iLenzuni et al] 
ll989l : IStasiriska & Szczerballl999l) . 

In our models, the spectrum of an AGN is approximated by 
a power law, i.e., L — CE~°', where E is the photon energy, L 
is the number of photons per unit time per unit of energy emitted 
by the source, and a is the spectral index. The constant C roughly 
indicates the radiation field intensity. It is usually defined in terms 
of the ionisation parameter, which is the ratio of the amount of 
hydrogen ionising photons incident on the inner border of the cloud 



U = 



Qt 



4:TTRf„nHC 



where uh is the numerical H density of the cloud, Rin is the in- 
ternal radius of the cloud, Qh is the number of hydrogen ionising 
photons emitted by the ionising source per unit of time, and c is 
the speed of light. U values between 10~^ an d 10~* are inferre d 
from observations of the NLR of AGNs (e.g. IOsterbrock|[l989l) . 
In this paper, we use models for average U values (0.1, 0.01 and 
0.001). The shape of the UV to X-ray SED can differ among AGNs 
( |Prietoetal]|2010l) . a range well represented by spectral indices 
from 1.0 to 1.5, which we use in our models. 

For comparison, we also include models with a starburst-like 
spectrum. For such models, we assume t hat the ionising spec- 
trum has the shape given by iRiffel et alj ( l2008bh . which repre- 
sents the underlying stellar population of the Starburst (SB) galaxy 
NGC 77 14. We scale the spectrum to study different number of ion- 
ising photons. T he models range from \og{QH) = 52 to 54, to 
cover SB values dLeitherer et al.|[l995l) . We study the same range 
of gas and dust parameters as for AGNs. 

Our models show that the bulk of the H2 near-IR emission in 
AGNs is produced in the warm region of the gas, so we use T > 
100 K as the stop criterion for the models. For SBs, there is a non- 
negligible contribution from gas with T < 100 K in some cases, so 
the stop criterion is T > 50 K. 



4.2 H2 Diagnostic Diagram: observations versus models 

Figure |6] shows the H2 line ratios resulting from our AGN models 
and from observations. Similarly, Fig. |7]compares our SBs models 
with observations. The observational data and the results of previ- 
ous models showed in Fig.[5]are also included in the plots of Figs. 
[6]and|7]for comparison. 

Our models with typical AGN parameters can reproduce most 
of the observed values (including the SFGs), while SB models only 
account for a fraction of the SFGs and BCGs and sub-sample of the 
AGNs with high 1-0 S(2)/l-0 S(0). Our models confirmed that the 
H2 excitation in AGNs is most likely produced in a gas irradiated 
by a spectrum rich in high energy photons, while for the SBs the H2 
excitation is could be due to a starburst-like spectra or by a power- 
law spectra as can be seen in Figs.|6]and|7] Models with a hard in- 
cident spectrum (i.e., with a large amount of high energy photons, 
or a lower a value) tend to have lower 2-1 S(l)/l-0 S(l) ratios than 
AGN models with soft spectra or SB models. The region around 
the coordinates (0.15,1.5) on these plots is densely populated with 
Seyfert galaxies and LINERs. For a =1.0, only models with 7ih < 
10" cm^'' can reproduce the observations. Models with a —1.5 
can reproduce the observations, with the exception of the dense 
(uh ~ 10^ cm~^) U =0.1 models. Alternatively, if the cloud is 
dense (nn ~ 10^ cm""'), the cloud is thiner than our models, i.e., 
has only gas at temperatures higher than 500K (see paragraph bel- 
low). The observations around the coordinates (0.12,1.5) can only 
be reproduced by models with nn ^ lO^cm^'*. It is worth men- 
tioning that the H2 line ratios remain essentially unchanged if we 
consider an average AGN ionizing spectrum like the one used by 
iKorista et al. or an equivalent power-law. For the latter we 

mean a SED that produces app roximately the same far-UV (5eV) 
to X-ray (2keV) flux ratio. The iKorista et al.l ( l 19971) SED, for ex- 
ample, is equivalent to a single power-law with a spectral index of 
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Figure 6. H2 diagnostic diagram: comparing AGN models and observations. Symbols represent observations, tlie polygons represent models found in the 
literature, and the dashed curve represents the ratios for an isothermal and uniform gas density distribution; they follow the same notation as in Fig. \5\ 
Calculated ratios are show in each panel as solid lines. Each panel shows AGN models for a different power law spectra, with the ionisation parameter U and 
the spectral index a as indicated. Each curve connects models for a dust grain size, ag = 10^^ (left curve) and 10^^/xm (right), where the gas density is 
varied. Dotted lines connect models of the same density for log n[f = 3, 4, and 5. These models assume that the cloud is at 1 pc from the central source. 
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Figure 7. Same as Fig.|6] but for starburst (SB) models. Calculated ratios are show in each panel as soUd lines, connecting models for different gas density. 
In each panel, the leftmost curve is for models with the dust grain size ag = 10^^ and the rightmost for 10~-^/im. Each panel shows SB models for different 
number of H ionising photons (Qh) ™d dust-to-gas ratio {M^/Mg). Qh is given in photons s^^. Dotted lines connect models of the same density for log 
nu = 3, 4, and 5. These models assume that the cloud is at 1 pc from the central source. 
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1.4. As an example, one of our models with a spectral index of 1.5 
gives the ratios 1-0 S(2)/ 1-0 S(0) = 1.252 and 2-1 S(l)/ 1-0 S(l) = 
0.523, while a model with the same parameters plus the Korista et 
al. (1997) SED gives very similar resuhs for 1-0 S(2)/ 1-0 S(0) = 
1.174 and 2-1 S(l)/ 1-0 S(l) = 0.566. 

The increase of the depth in the cloud also affects significantly 
the line ratios as can be seen in Fig. [8] In general, for increasing 
depth both ratios decrease towards lower values, initially following 
the LTE high temperature solutions (dashed curve) and then fol- 
lowing a path that depends on object properties. It can clearly be 
observed in this figure that the ratios are significantly affected if 
the models run only for T ^ 1000 K (marked with red stars), thus, 
affecting model predictions. 

The diagrams show that several factors, besides the radiation 
source, compete for the determination of H2 line ratios. Gas and 
dust density, grain size, and the shape and luminosity of the central 
source spectra are the most important. For example. Figure |6]shows 
that AGN models with high gas density tend to have low 1-0 S(2)/l- 
S(0) and 2-1 S(l)/l-0 S(l) ratios, reflecting an emission spect rum 
domin ated by collisions, as discussed by Sternberg & Dalg arnd 
( Il989l) . Values obtained bv lKwarj j 19771) for shock regions are sim- 
ilar to our densest AGN models. For SBs, high density models also 
tend to have low values of 1-0 S(2)/l-0 S(0), but the exact behaviour 
of 2-1 S(l)/l-0 S(l) with density also depends on other properties, 
as discussed below. 

Since we assume a constant dust-to-gas ratio (and not the den- 
sity of grains), increasing the gas density also increases propor- 
tionally the amount of dust. As seen in Fig. [T] increasing just the 
dust-to-gas ratio (keeping nn constant) diminishes the 1-0 S(2)/l- 
S(0) ratio, but not necessarily 2-1 S(l)/l-0 S(l). We show only 
SB models with different dust densities, but the effect of increasing 
dust density on AGN models is similar to that on SB models. 

As previously mentioned, dust grain size also plays an im- 
portant role in H2 line ratios. Models with smaller dust grains 
have smaller 2-1 S(l)/l-0 S(l) ratios. Assuming a constant dust- 
to-gas ratio, the optical depth due to dust and the H2 forma- 
tion rate are in versely proportional to grain size, as shown in 
lAleman & Grue nwald (2004). For these reasons, more H2 is pro- 
duced in models with smaller grains. 

4.3 H2 excitation temperatures and masses 

As stated in iRodriguez-Ardila et al. 1 ( l2005l) . an alternative way to 
determine the mechanism driving the molecular gas excitation is 
to derive the rotational and vibrational temperatures for the H2. 
A gas dominated by thermal excitation has similar rotational and 
vibrational temperatures (as expected for a gas in LTE) while 
fluorescent excitation is characterized by a high vibrational tem- 
perature and a low rotationa l temperature. For more details see 
iRodriguez-Ardila et al.l ( I2OO5I and references therein). 

The values of Tvib and Trot can be calcu- 
lated using the fluxes of the observed H2 lines to- 
gether with the expressions 
Tvib = 5600/In(1.355 x 



Table 5. Extinction coefficient. Moleculai' gas mass, vibrational and rota- 
tional temperatures. 



-1113/ln(0.323 x 7i_ 



l lReunanen et al.l I2OO2I) : 

-^i-os(i)/-^2-is(i)) and 



US(2)/-'l-0S(0); 



The derived 



values ar e presented in Table [5] It can be observed from that table 
and from iRodrfguez-Ardila et aL ] 1 I2OO5I) Table 4, that AGNs and 
LINERs tend to have similar values of Tvib and Trot, while in 
the case of SFGs (NGC 23, NGC 520, NGC 34, NGC 1614 and 
NGC 7714, the latter 3 are from Rodrigues-Ardila sample) the 
Tvib tends to be higher than Tot . 

In addition to the excitation temperature, the fluxes of Ta- 



Source 


Cext 


H2 Mass 


Tvib 


Trot 






MO 


K 


K 


NGC 23 


- 


- 


2521 ±469 


1067±601 


NGC 520 


- 


- 


4115±279 


1832±159 


NGC 660 


0.69 


161±5 


2704±240 


2108±208 


NGC 1055 


- 


- 


- 


- 


NGC 1134 










NGC 1204 


0.59 


1991±57 


2912±256 


2776±270 


NGC 1222 


0.21 


67±7 


< 2856 


2373±399 


NGC 1266 


1.35 


7189±132 


2232±154 


1631±118 


UGC 2982 










NGC 1797 


0.54 


1606±56 


4293±271 


1209±307 


NGC 6814 


0.00 


56±19 


3607±381 


1165±343 


NGC 6835 


0.98 


201±6 






UGC 12150 


0.43 


3227±97 


2444±363 


960±630 


NGC 7465 


0.47 


323±34 


2362±227 


1378±222 


NGC 7591 


0.51 


2638±82 


3478±361 


1069±174 


NGC 7678 


0.28 


170±17 


4225±806 





ble m also allow us to compute the mass of hot H2 emitting the 
l-OS(l) 2.121 fim line by means of the e quation: — 5.0875 x 
10^^T»^/i_os(i) ( iReunanen et alj2002h . 

The mass of this hot H2 for our galaxy sample is listed in Ta- 
ble[5] They were calculated ass uming T=2000 K, a transition prob- 
ability As(i)=3.47x 10"'' s"^ ( lTumeretai1ll977h . the population 
fracti on in the i/=l, J=3 level /^^i, j=3=0.0122 jScoville et al.l 
Il982h and the intrinsic flux of H2 2.121^m, /i_os(i)- The extinc- 
tion coefficient, Cext, is calculated by assuming an intrinsic value 
of 5.88 for the flux ratio Pa/3/Br7 1 Osterbrock 1989 , , case B). A s 
can be observed in Table |5] and in Rodriguez-Ardila et al.l (|2005|) . 
the mass of the hot H2 is very similar for all activity types. The 
fraction of molecular mass present in the nuclear region and emit- 
ting in the NIR is a very small fraction of the warm molecular 
mass expected to be present in the galaxy centre (up to lO^^Af©, 
lYoung et al.lll96lL for example), and is not related to activity type. 

Overall, the diagnostic diagrams presented in Figs. [6] and |7] 
allow us to distinguish fairly well between the H2 excitation mech- 
anisms (collisional or UV pumping). However, it is not straightfor- 
ward to distinguish between different classes of sources (e.g. AGN 
or Starburst), as discussed in the previous subsection. Alternatively 
mechanisms may dominate the H2 excitation in different regions 
in a same object. As the emission lines are integrated along the 
line of sight, we observe a combination of such processes. As can 
be seen in Figs. |6| and |7] both collisional and UV pumping exci- 
tations may play an important role not only in AGNs but also in 
SFG galaxies. This is further supported by the similarity between 
the vibrational and rotational temperatures of H2 in some objects, 
and the tendency of Tvib to be higher than Tot in others (see Ta- 
blejS)- For example, in SFGs, the thermal processes may dominate 
the excitation of the gas irradiated by supernova remnants (SNRs) 
or very hot stars. In this case, the H2 emission is due to X-ray exci- 
tation, whi ch would favour vibrational tr a nsitions over r otational 
transitions ( iRodrfguez-Ardila et al. I I2OO5I) . iMouril ( I1994 showed 
that SNRs lies very close to the thermal curve in these diagnostic 
diagrams (see their Fig. 1). In the case of UV excitation by stars, 
the effective optical depth of the nebula, a non-thermal signature, 
will be determined if the ratios reflect the domination of thermal or 
non-thermal process (Sternberg & Dalgarno 1989). 
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Table 4. Fluxes of atomic and molecular lines for our sample, in units of 10 ergs cm ^ s ^, measured in the Sample. 



Source 


[Fen] 


Pa/3 


[Fell] 


H2 


H2 


H2 


Ha 


H2 


Br7 




1.2570/im 


1.2820/im 


1.6444/im 


1.9570/^m 


2.0332/im 


2.1213/im 


2.2230/im 


2.2470/^m 


2.1650/im 


NGC 23 


10.50±0.87 


— 


14.20±3.46 


20.20±4.60 


4.94±0.57 


10.00±1.26 


4.53±2.10 (t) 


1.47±0.37 


0.00±0.00 


NGC 520 


— 


— 


3.87±0.14 


2.41±0.39 


1.13±0.10 


2.40±0.18 


0.67±0.06 


0.83±0.12 


6.67±0.19 


NGC 660 


13.90±0.65 


29.0±0.60 


15.70±0.90 


8.87±0.80 


3.67±0.56 


6.91±0.77 


<2.01 


1.18±0.12 (t) 


16.60±0.71 


NGC 1055 


— 


— 


— 


— 


O.OOiO.OO 


— 


— 


— 


— 


NGC 1134 










O.OOiO.OO 










NGC 1204 


5.16±0.2 


12.10±0.20 


6.51±0.36 


4.17±0.51 


1.81±0.34 


3.61±0.25 


0.87±0.09 


0.71±0.10 


5.88±0.27 


NGC 1222 


6.32±0.43 


28.20±0.37 


4.65±0.18 


1.95±0.34 


0.95±0.14 


0.84±0.12 


0.49±0.14 


<0.16 


6.98±0.073 


NGC 1266 


3.19±0.51 


0.75±0.08 


2.90±0.38 


15.10±0.47 


5.10 ±0.44 


13.70±0.25 


3.26±0.12 


1.51±0.14 


1.40±0.33 


UGC 2982 








1.79±0.48 


0.61 ±0.08 (t) 


0.60±0.02 (t) 






0.80±0.036 


NGC 1797 


5.04±0.25 


12.10±0.26 


4.65±0.47 


4.10±0.51 


1.48±0.30 


3.21±0.28 


1.20±0.17 


1.18±0.16 


5.33±0.09 


NGC 6814 


4.64±0.52 


3.53±0.57 


5.53±0.48 


3.36±0.54 


1.25±0.20 


2.14±0.20 


1.05±0.23 


0.61±0.13 


0.52±0.32 


NGC 6835 


2.06±0.23 


4.86±0.19 


3.24±0.08 


3.33±0.44 


0.85±0.05 


1.34±0.21 


0.65±0.18(t) 




4.66±0.25 


UGC 12150 


4.73±0.33 


7.86±0.30 


3.56±0.18 


5.25±0.05 


1.65±0.07 


4.00±0.13 


1.70±0.85 (t) 


0.55±0.12(t) 


2.88±0.13 


NGC 7465 


11.8±0.66 


9.57±2.61 


9.20±0.28 


4.33±0.16 


2.36±0.39 


3.92±0.27 


1.71±0.10 


0.49±0.06 


3.75±0.69 


NGC 7591 


6.80±0.21 


9.74±0.21 


6.05±0.60 


8.17±0.26 


2.46±0.20 


4.80±0.36 


2.25±0.25 


1.30±0.26 


4.07±0.052 


NGC 7678 


3.28±0.49 


9.20±0.51 


3.09±0.21 


1.45±0.10 


2.03±0.20 


0.88±0.14 


>0.71 


0.32±0.14 


2.56±0.15 



(t) Affected by telluric absorption. 




the e nd of Seyfert-lilce nuclei and the begi nning of LINER-like ob- 
iects. lRodriguez-Ardila et aL I( l2004l2005h confirm such diagram as 
a suitable means of separating emission-line objects by their degree 
of activity in the NIR, and propose that AGNs are characterized by 
H2 2.121/im/Br7 and [Fell] 1.257^m/Pa/3 flux ratios between 0.6 
and 2. Starburst/H II galaxies display line ratios <0.6, while LIN- 
ERs are characterized by values higher than 2 for both ratios. How- 
ever, the lack of an adequate number of objects preclude definitive 
conclusions in those studies. 

Here, with a larger sample, and most importantly, a more ade- 
quate number of objects of each activity class, we may provide fur- 
ther support to the trend already observed. Figure |9]presents the up- 
dated version of H2 2.121/im/Br7 x [Fell] 1.257/im/Pa/3 diagnos- 
tic diagram. The sample analyzed here is composed by the objects 
listed in Tabled] plus some taken fr om the liter ature: AGNs from 
[Rodrig uez- Ardila et al. I ( |2004 l2005l) . SFG from bale et al.l ( l2004h 



Figure 8. Behavior of the ratio 1-0 S(2)/l-0 S(0) vs 2-1 S(l)/l-0 S(l) with 
the distance from the inner border of the cloud, for several AGN models. 
Each solid curve represents a model, with distance increasing from right to 
left. The red stars and black dots over these curves indicate the stop criterion 
of T = 1000 K and 100 K, respectively. Dots are observations and the dashed 
curve represents the ratios for an isothermal and uniform gas distribution as 
shown in Figs.[6]to[2] 



5 NUCLEAR ACTIVITY FROM NEAR-INFRARED 
EMISSION-LINES 

Since H2 and [Fell] lines are common features in a wide vari- 
ety of sources, it is interesting to investigate how these emission- 
lines compare in objects with different degrees of nuclear ac- 
tivity. A diagram involving the line rat ios H2 2.121 / im/Br 7 and 
[Fell] 1.257^m/Pa/3 was proposed by iLarkin et al. ] ( Il998h m a 
study of LINERs and other emission-line objects, including a few 
Seyferts. They report a strong linear correlation in the log-log plot 
of [Fe Il]/Pa/3 versus H2/Br7, with SB galaxies displaying the lower 
values, Seyferts with intermediate ones and LINERs with the high- 
est ratios. They suggest that [FeIl]/Pa/3 ~1 and H2/Br7 ~3 mark 



and iL arkin et al. ( 1998|) . LINERs, supe rnova remnants (SNRs ) 
from lLarkin et al I ( Il998l) . and BCGs from llzotov & ThuanlJioTll) . 
As can be seen in Figure|9] LINERs and AGNs are not clearly sep- 
arated by such a diagram. However, our sample of LINERs is com- 
posed by Luminous Infrared Galaxies (LIRGs), and the classifica- 
tion in Table [T] is based on opti cal emission-line diagrams, w hich 
are affected by reddening (see iRodriguez-Ardila et al.ll20oH for 
example). Since LIRGs harbour large quantities of dust, we sug- 
gest that the LINERs lying within the region of AGNs may hide a 
Seyfert nucleus. 

With our sample comprising ~ 65 objects, we not only con- 
firm the diagram (Fig. |9) as a discriminator of emission-line ob- 
jects by their degree of activity, but we also provide improved lim- 
its to the line ratios for each activity type. The limits are as fol- 
lows: (i) SFGs [Fell] 1.257^m/Pa^ < 0.6 and H2 2.121^tm/Br7 < 
0.4, (ii) AGNs 0.6 < [Fell] 1.257^m/Pa/3 < 2 and 0.4 < H2 
2.121/.im/Br7 < 6, and (iii) LINERs [Fell] 1.257^m/Pa^ > 2 and 
H2 2.121/im/Br7 > 6, 

As discussed in [Rodriguez- Ardila et al.l ( l2005h on a pure ob- 
servational basis, the interpretation of Figure |9] is not straightfor- 
ward. For example, the correlation between H2/Br7 and [Fe Il]/Pa/? 
may originate from the fact that H2 and [Fe II] are mainly excited 
by the same mechanism, i.e.. X-ray. In fact, it was shown very re- 
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cently by iDors et al] ( |2012|) . using photodissociation models, that 
X-ray heating from the AGN is an important mechanism in the ex- 
citation of these NIR emission-lines in AGNs (see Sec|4l). They 
also show that a decrement in the X-ray content of the contin- 
uum source translates into a line weakening, and their models are 
no longer compatible with observations. However, X-ray heating 
is a plausible exp l anatio n for AGNs, but may not work for LIN- 
ERs. lLarkin et all ( Il998l) . for example, argue that hard X-ray heat- 
ing from a power-law is a plausible mechanism to explain LINERs 
with low values of the [Fe Il]/Pa/3 and H2 2.121/^m/Br7 ratios. But, 
it would not explain the high values (up to ^-^7) observed in some 
objects. Nevertheless, the high [FeIl]/Pa/3 line ratios may be ex- 
plained by the dependence of [FeIl]/Pa/3 on the Fe/O abundance 
and on the dependence on the ionisation para meter o f the molecu- 
lar gas, leading to high values for both ratios dOors e t al. 2012, see 
their Figure 1) 

On the other hand, it has been recently proposed that it is un- 
likely that LINER emission lines are powe red by photoionisation 
by the low-lum inosity A GN they harbour ( Eracleous et al.ll2010l ; 
ICid FemandeTe t al. 201(j: lYan & Blantonll2012h . This would rule 
out X-ray heating as the main mechanism driving the excitation of 
the NIR emission-lines in LINERs. We investigated the role played 
by the hard X-ray (2-10 KeV, F^;) radiation in [Fe II] and H2 emis- 
sion lines with diagrams involving F^fBi^ vs H2/Br7 ratios and 
Fx/Psl/B vs [FeIl]/Pa/3 (Fig.llOt. The X-ray fluxes were taken from 
literature and are listed in Table |6] A mild correlation is observed 
for Fx/Bi"/ vs H2/Br7 (Fig. 1 10b). indicating that the X-ray heating 
is an important mechanism in the molecular gas excitation. For a 
more objective assessment of the X-ray heating in our data, we in- 
cluded model predictions to Fig. llOb . It is clear that the explored 
parameter space overlaps the observed line ratios. In addition, the 
AGNs are better explained by a=l and (7=0.1. It is also interesting 
that Sy 1 s require higher hydrogen gas densities than the other ac- 
tivity types, indicating that the molecular gas in Sy Is is denser. No 
correlation is observed for Fx/Pafi, and, [Fe Il]/Pa/3 and as for the 
diagram invloving H2/Br7 no separation is observed between LIN- 
ERs and Seyferts (Fig. [TOt . This result suggests that X-ray heat- 
ing alone is not respons i ble fo r the observed sequence in Fig. |9] 
As shown bv lDorsetaLl ( l20I2l) . high values of [Fe Il]/Pa/? line ra- 
tios can be explained by the dependence on the Fe/O abundance. 
Thus, the observed sequence in Fig.|9]may be a consequence of the 
amount of X-ray photons together with abundance effect (i.e. in the 
case of LINERs powered by star-formation, SN explosions would 
eject enriched material to the interstellar medium). 

A detailed interpretation of the variation of these line ratios, 
from pure-starburst driven emission to LINER nuclei, is non-trivial. 
A possible scenario is that the different line ratios reflect different 
excitation mechanisms, in the sense that they are able to describe 
a transition from an ionizing radiation p owered by star forma tion 
to pure shock excitation driven by SNRs jMazzalav et al.l2012h . In 
such a case, the emission-line ratios observed in our LINER sample 
would be a combination of processes, from X-ray heating to star- 
formation/evolution. Since most LINERS have strong star forma- 
tion (Tab.[T]l, the stars could provide the required electromagnetic or 
mechanical power to produce the observed line strengths. However, 
to clearly investigate this issue, a 3D simulation of the emitting 
gas would be needed (i.e. AGN in the center -1- star-formation in a 
volume). It is also worth mentioning that, in the case of star forma- 
tion, the underlying stellar population would significantly affect the 
atomi c hydrogen em ission-lines, thus affecting the observed line 
ratios jLarkin et al.iri 998: Riffel et al. 2008bj. 

Another interesting point is the fact that Br7 is produced 



Table 6. Hard X-Ray Fluxes, in units of 10 
from literature. 



taken 



Source 


Fx(2-10 KeV) Reference 


Sy 1 


Mrk 334 


80 


Lutz et al. (2004) 


NGC 7469 


290 


Lutz et al. (2004) 


NGC 3227 


75 


Lutz et al. (2004) 


NGC 4151 


900 


Lutz et al. (2004) 


Mrk 766 


300 


Lutz et al. (2004) 


NGC 4748 


30 


Landi et al. (2010) 


NGC 5548 


430 


Lutz et al. (2004) 


NGC 6814 


300 


Turner & Pounds (2010) 


Sy 2 


ESO428-G014 


3.8 


Lutz et al. (2004) 


NGC 591 


2 


Guainazzi et al. (2005) 


Mrk 573 


2.8 


Guainazzi et al. (2005) 


Mrk 1066 


3.6 


Guainazzi et al. (2005) 


NGC 21 10 


430 


Nandra et al. (2007) 


NGC 7674 


5 


Lutz et al. (2004) 


NGC 5929 


79 


Lutz et al. (2004) 


Mrk 1210 


250 


Lutz et al. (2004) 


NGC 5728 


13.3 


Shu et al. (2007) 


Mrk 993 


6.1 


guainazzi et al. (2005) 


NGC 5953 


0.06 


Guainazzi et al. (2005) 


NGC 1144 


1100 


Lutz et al. (2004) 



LINER 



NGC 5194 


4.8 


QiPDi et al. (2006) 


NGC 7743 


0.44 


Gonzalez-Martin et al. (2009) 


NGC 660 


0.13 


Filho et al. (2004) 


NGC 7465 


41.5 


Guainazzi et al. (2005) 


NGC 3998 


84.2 


Uedaetal. (2001) 


NGC 7479 


3.9 


Ueda et al. (2005) 


NGC 4736 


19 


Uedaetal. (2001) 


SFG 


NGC 34 


39 


Lutz et al. (2004) 


NGC 7714 


2.6 


Ueda et al. (2005) 


NGC 1614 


3.6 


Ueda et al. (2001) 


NGC 3310 


17.4 


Uedaetaf (2001) 


BCDg 


Mrk 930 


0.161 


Gonzalez et al. (2009) 



mostly inside the hydrogen ionised region, while the H2 emission 
origin ates in a warm semi-ionised region dAleman & Gruenwaldl 
I2OI If) around the hydrogen recombination and in the photodissoci- 
ation region (PDR), where H is mostly neutral. SFGs typically have 
lower 1-0 S(l)/Br7 ratios than AGNs. This may indicate that the re- 
gion that produces the H2 emission in AGNs is larger than in SBs, 
if we consider that AGN spectra are harder than SFG's. Harder ion- 
isation spectra may produce a more extended semi-ionised region 
(Aleman & Gmenwald 2004, 2007), which favours the formation 
and emission of H2. In our models, the ratio H2 1-0 S(l)/Br7 is 
within the range 10"^ <H2 1-0 S(I)/Br7 < 10 (see SeclU, with 
the higher ratios (> 1) occurring in models with a —l.O and the 
smallest values (< lO"'^) occurring for a =1.5 (in fact, the emis- 
sion of these last objects are below the detection limit of our obser- 
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vations). Furthermore, denser and/or dustier models typically have 
higher ratios. A more detailed discussion on these models will be 
made in a forthcoming work (Aleman et al., in preparation). 



6 CONCLUDING REMARKS 

We present NIR spectroscopy (0.8/im — 2.4/im) of seven SFGs 
and nine LINERs to discuss the distribution and excitation of the 
emitting H2 and [Fell] gas in emission-line galaxies. These data 
were anal ysed together with similar AGN data of two previous pub- 
lications iRodn'guez-Ardila et ani2004l2005h . The whole data set 
constitutes the most complete and homogeneous sample of such 
kind of objects observed in the NIR to date. Our main conclusions 
are: 

• H2 is common within the inner few hundred parsecs of 
emission-line galaxies, regardless of type. However, the molecu- 
lar gas follows different kinematics than that of the ionised gas, 
suggesting that the two emissions are not co-spatial. For example, 
it is clear that the H2 lines are narrower than the forbidden lines, 
especially for [S III]. This could be interpreted as if the latter is 
broadened by the gravitational influence of the SMBH in the AGNs. 
Furthermore, LINERs tend to have [S III] significantly broader than 
[Fe II] and H2 compared to SFGs. The fact that H2 FWHM is unre- 
solved in almost all objects implies that the molecular gas is proba- 
bly not gravitationally bound to the SMBH but to the gravitational 
potential of the galaxy following a different kinematics than that of 
the classical NLR gas. In addition, IFU data and rotation curves de- 
rived for H2 and published by other authors for objects in common 
with our sample support this hypothesis and point to a scenario 
where H2 may be arranged in a disc-like structure on the galaxy 
plane. 

• We have computed new photoionisation models for clouds 
ionised by AGN and stars, which provided fluxes of the H2 lines. 
We use a diagnostic diagram with H2 IR line ratios to study the ex- 
citation mechanisms of H2. The thermal excitation plays an impor- 
tant role not only in AGNs, but also in SFG galaxies. This hypoth- 
esis is further supported by the similarity between the vibrational 
and rotational temperatures of H2 in some objects, and the tendency 
of Tvib to be higher than Trot in others. In SFGs, the importance 
of the thermal excitation may be associated with the presence of 
SNRs close to the region emitting H2 lines. 

• With our extended sample, of 65 objects, we have confirmed 
that the diagram involving the line ratios H2 2.121/im/Br7 and 
[Fe II] 1 .257/im/Pa/3 is an efficient tool for separating emission-line 
objects according to their dominant degree of activity and should 
be useful for classifying objects with hidden AGNs or highly red- 
dened objects. We also provide improved limits to the line ra- 
tio intervals for each activity type, as follows: (i) SFGs [Fell] 
1.257^m/Pa/? < 0.6 and H2 2.121^m/Br7 < 0.4, (ii) AGNs 0.6 
< [Fell] 1.257^m/Pa/? < 2 and 0.4 < H2 2.121/im/Br7 < 6, and 
(iii) LINERs [Fell] 1.257^m/Pa/J > 2 and H2 2.I21/im/Br7 > 6, 

• A positive correlation is found in the diagnostic diagram in- 
volving the above ratios, interpreted as being a more likely an ap- 
parent correlation, since different mechanisms drive the production 
of either [Fe II] or H2 according to the level of nuclear activity. It 
is not discarded that a combination of X-ray heating from the cen- 
tral source (in AGNs) plus shock heating (radio jet and NLR gas 
interaction) and/or SNRs can simultaneously drive the emission of 
[Fe II], supporting this correlation in LINERs and SFGs. 
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Figure 9. Diagnostic diagram: boxes are AGNs, Sy 1 in black (open) and Sy 2 in red (filled), asterisks are LINERs (green), filed circles are SFG (blue), triangles 
(magenta) are SNRs and open circles are blue compact dwaif galaxies (violet). The solid line is [FeII]/Pa/3 = 0.749(±0.072)H2/Br7 — 0.207 it 0.046 with 
correlation coefficient = 0.80 and r^=0.64. 
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